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We study scaling behavior of the geometric tensor Xa,p{Xi, X2) and the fidelity susceptibility (xf) 
in the vicinity of a quantum multicritical point (MCP) using the example of a transverse XY model. 
We show that the behavior of the geometric tensor (and thus of xf) is drastically different from that 
seen near a critical point. In particular, we find that is highly non-monotonic function of A along the 
generic direction Ai ~ A2 = A when the system size L is bounded between the shorter and longer 
correlation lengths characterizing the MCP: l/|A|' yi <C L <C 1/|A|" 2 , where v\ < V2 are the two 
correlation length exponents characterizing the system. We find that the scaling of the maxima of 
the components of Xafi is associated with emergence of quasi-critical points at A ~ 1/L 1 ^ 1 , related to 
the proximity to the critical line of finite momentum anisotropic transition. This scaling is different 
from that in the thermodynamic limit L 3> 1/|A|" 2 , which is determined by the conventional critical 
exponents. We use our results to calculate the defect density following a rapid quench starting from 
the MCP and show that it exerts a step-like behavior for small quench amplitudes. Study of heat 
density and diagonal entropy density also show signatures of quasi-critical points. 
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I. INTRODUCTION AND MAIN RESULTS 

Fidelity and fidelity susceptibility^— (xf) are infor- 
mation theoretic measures of continuous quantum phase 
transitions (QPT) . It is well known that a QPT&i, oc- 
curring at absolute zero temperature, indicates a change 
in symmetry of the ground state of a many body quan- 
tum system. Fidelity measures the overlap between two 
neighboring ground states in the parameter space of a 
quantum Hamiltonian and shows a dip at the quan- 
tum critical point (QCP) where the overlap between two 
neighboring ground states is minimum. At the same 
time, the fidelity susceptibility, which is a quantitative 
measure of the rate of change of the ground state of a 
system under an infinitesimal variation of a parameter of 
the Hamiltonian, usually attains the maximum value and 
often diverges with the system size in the thermodynamic 
limit. The scaling behavior of xf with the system size 
at the QCP and with respect to the deviation from the 
QCP is given in terms of the associated quantum critical 
exponents. 

Following the Kibble Zurek (KZ) predictions^— on the 
scaling of the defect density as a function of rate of change 
of the Hamiltonian following a slow quantum quench (see 
Refs. [HUH! for a review), a series of works have been 
reported addressing the scaling behavior of the defect 
density for quench across a QCP, gapless phases, gap- 
less lines, etoi£. Of a particular interest is generalization 
of these results to quenching across a multicritical point 
(MCP)i^~— . Such points (both for classical and quan- 
tum phase transitions) are characterized by two or more 
orthogonal directions with different correlation lengths^. 
They generically appear either as intersection of two crit- 
ical lines or as a terminating point of one critical line on 
the other critical line. The second scenario is realized 
for the transverse field XY model, which is the primary 
example for this paper (see Fig. [TJ . In this model one 



terminating critical direction (horizontal line in Fig. [T]) 
corresponds to a finite momentum anisotropic transition 
with the critical momentum vanishing at the MCP. In 
finite size systems the momentum is quantized in units 
of 1/L. In turn this quantization results in emergence 
of quasi-critical points close to the MCP where the en- 
ergy gap has a local minimum. As we will show these 
quasi-critical points dominate the scaling of the fidelity 
susceptibility xf along generic directions and in turn de- 
termine the scaling of defect density for slow as well as 
rapid quenches across a QCP as suggested in Refsi22r— . 

The scaling behavior of xf at a QCP is well 
established 2 ^—. For a Hamiltonian given by H = Hq + 
XHj = Hq + XV where Hq is the Hamiltonian describing 
the QCP, V = c?a^| A _ = J2r v ( r ) i s the perturbation 
not commuting with Hq and A measures the deviation 
from the QCP. For the translationally invariant systems 
it can be shown that the fidelity susceptibility is related 
to the connected imaginary time correlation function of 
the perturbation^: 

X f(X)=J dr J drr(v(r,r)v(0,0)) c . (1) 

From this representation it is clear that the scaling di- 
mension of xf is dim[xF] = 2A,j — 2z — d. If this scaling 
dimension is negative then the fidelity susceptibility di- 
verges at the QCP as xf(A = 0) ~ L d+2z ~ 2A ", where 
L is the system size. If the scaling dimension is positive 
then the fidelity susceptibility is still singular but this 
singularity is a sublcading correction on top of a nonuni- 
versal constant (see e.g. Ref. [H[). In case where V 
is a marginal or relevant perturbation (so Xv(r) scales 
as the energy density) we have additional simplification 
coming from A.„ = d + z — \ jv so that at the critical 
point^^ 

XF ~ L 2 /»~ d . (2) 
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At small but finite A this scaling is substituted by 

Xf ~ |Ar d - 2 (3) 

and the two asymptotes match when the correlation 
length £ ~ becomes of the order of the system 

size. These asymptotics are dominant for dv < 2 and 
subleading for dv > 2 (at dv = 2 there are additional 
logarithmic singularities 2 ^) . 

It is natural to expect that the scaling of the fidelity 
susceptibility can be further extended to understand uni- 
versal properties of MCPs as well as quantum dynamics 
through the MCPs. For example, in Refs. (22ll23j it was 
shown that the scaling of the generated quasiparticlcs 
and excess energy (heat) produced by slow quenches can 
be understood through the scaling behavior of general- 
ized adiabatic fidelity susceptibilities. To the best of our 
knowledge such analysis was never performed before. As 
we mentioned earlier because MCPs are characterized by 
at least two orthogonal directions (principal axes) char- 
acterized by different scaling exponents, the natural ob- 
ject to analyze near the MCP is the so called geometric 
tensor— defined as 



X«4» = y>|0A a |n><n|3xJO>, 



n#0 



(4) 



where the sum is taken over all excited states. One can 
also define the fidelity susceptibility, characterizing the 
drop of the overlap of the ground state wave functions for 
an infinitesimal change of the tuning parameter 8 A. It 
is clear that now \f is described by both the magnitude 
and the direction: 



a, f) 



Xap rip, 



(5) 



where n a is the projection of the vector <5A to the direc- 
tion a. If we are interested in the fidelity susceptibility 
along a generic direction such that all components n a are 
of the same order then the most divergent component of 
the geometric tensor defines the scaling of xf- It fol- 
lows from the definition of the geometric tensor that its 
diagonal components describe the fidelity susceptibility 
along the principal axes. For the purposes of this work 
we will focus on the situation of two orthogonal direc- 
tions a, f3 — 1,2. It is straightforward to generalize our 
results to the situation of multi-component MCPs. 

It is natural to expect that both the geometric ten- 
sor and fidelity susceptibility will have singular behavior 
at the MCP. In this work we focus on how this singu- 
larity develops as we approach the MCP. We find that 
this approach to the MCP is highly nontrivial charac- 
terized by non-monotonic dependence of the components 
of the geometric tensor on the coupling A. In particu- 
lar, in the thermodynamic limit L 3> £i(A),£2(A), where 

6.2(A) 

sor can be expressed as: 



|Ai,2| iyi ' 2 , we find the that the geometric ten- 



where d is the dimensionality, v\ and v^ are the corre- 
lation length exponents corresponding to the couplings 
Ai and A2 respectively, and F a .p is a scaling function. 
For approaching the MCP from a generic direction Ai = 
Asin6* and A2 = Acos# with A — > and 6 7^ 0,7r/2, .. 
fixed, the scaling of the geometric tensor is affected by 
the small argument asymptotic of the scaling function 
Fa,p{x) at x — > (we assume that v\ < v^). While for ap- 
proach of the critical point along the line IA1I" 1 cx IA2I" 2 
the argument of the scaling function is constant of the 
order of unity and the scaling of the geometric tensor 
is determined purely by the prefactor in Eq. ©. We 
point that some asymptotics of the scaling function in 
Eq. ([6]) can be constrained by critical exponents char- 
acterizing second order phase transitions away from the 
MCP. Thus if the A2 = describes a line of critical points 
then the asymptotics of X22 at IA1I" 1 IA2I" 2 should re- 
duce to the scaling behavior of fidelity susceptibility of 
the corresponding phase transition (see Eq. (j3])). Clearly 
this requirement determines the asymptotics of i*22(x) 
at small x. One can find similar constraints on other 
components of F a ^{x) at large and small x. In the next 
section we analyze the scaling functions and discuss their 
asymptotics for a particular one-dimensional transverse 
field XY-model. 

In finite size systems the scaling (|6|) is replaced by an- 
other expression: 



Xa,/j(A,L) ~ l}l"<*+^-*F ajP {L\\ x \ v \L\\ 2 \^). 



(7) 



(6) 



If L > 1/ 1 Ai I Ul , 1/ 1 A 2 1 " 2 then the scaling above should re- 
duce to Eq. ©. Exactly at the MCP (A = 0) the geomet- 
ric tensor should be independent on A which leaves only 
two possibilities F a ^(0, 0) = const > or F a ,p(0, 0) = 0. 
Thus the nonzero components of the geometric tensor 
should have the scaling, which is a straightforward gen- 
eralization of Eq. ([2]). As it is clear from Eq. ([7]) there 
is a third scaling regime where one argument of the 
scaling function F is small and the other is large. For 
generic direction Ai ~ A2 = A this regime happens 
when the system size is bounded between the shorter 
and longer correlation lengths characterizing the system: 
l/|A| yi < L < 1/1 A" 2 1 with a new scaling emerging 
when L w £1 ~ l/IA^ 1 . For the model we analyze we 
find that all components of the geometric tensor have 
additional singularity at this (quasi-critical) point lead- 
ing to the local maximum of Xa,/3 with stronger scaling 
with L. We also find that in this intermediate regime 
1 / 1 A I ^ x < L < 1 / 1 A I ^ 2 all components of the geometric 
tensor have strongly oscillatory dependence on the mag- 
nitude of A (or equivalently on the system size). These os- 
cillations, however, maybe an artifact of our model where 
MCP corresponds to the ending point of a nonzero mo- 
mentum anisotropic transition. And finally at L ~ £2 (A) 
the two scalings © and (JT]) smoothly connect to each 
other. 

The oscillations of the components of the geometric 
tensor and the fidelity susceptibility near the MCP in 
turn translate into a step-like behavior in the number of 
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defects produced during a sudden quench of the magni- 
tude |A| starting from the critical point. The steps are 
more pronounced in relatively small systems, while in 
the limit of large system sizes we get a scaling behavior 
n cx ~ |A|' yid . This result is similar to that observed in 
the rapid quench from a critical poinli^ 2 .. 



II. GEOMETRIC TENSOR FOR THE 
TRANSVERSE FIELD XY-MODEL 

Let us begin with the spin- 1/2 one-dimensional trans- 
verse XY Hamiltonia n 10 ! 26 " — 



MC . 



J + J. 
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y 
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(8) 



where cr's are the Pauli spin matrices, L is the system 
size, J x , J y are the coupling constants along x and y 
directions respectively, and h denotes the transverse field. 
Here we consider even boundary condition {c p L+1 = <r v L , 
where p = x, y, z) in the odd fermion number sector, i.e., 

P = | ^1 — n^=-L °j) = 1- As wc will explain later 
choice of boundary condition does not affect the results in 
the thermodynamic limit. The boundary conditions will 
somewhat modify the finite size scalings (see below) but 
the exponents will remain unaffected. It is well known 
that the Hamiltonian (8) can be mapped to a collection 
of non-interacting quasipartices, labeled by momentum 
modes k (fc = 0, n/L, 2ir/L, 27r(L — 1)/L), with an 
excitation energy of the for m 30 i 31 



A fc = y/(h+ Jcosk) 2 +7 2 sinfc 2 , 



(9) 



where J = (J x + J y ) and 7 = (J x — J y ). Clearly the 
spectrum is gapless along the line 7 = 6 provided that 
—l<h/J<l and along the two lines h = ±J. At 
these boundaries the system undergoes continuous quan- 
tum phase transitions. The corresponding phase diagram 
is plotted in Fig.[TJ The points where the two critical lines 
meet arc multicritical. We will focus on the MCP where 
J = h and 7 = (denoted by MCI on Fig. [TJ. Then 
the vanishing gap corresponds to the momentum k = n 
therefore it is convenient to shift momentum k — > tt — k. 
Then the spectrum (9) linearized near the MCP simpli- 
fies to 



A fe a + Jfc 2 /2) 2 + 7 2 fc 2 (10) 

The ground state of the Hamiltonian ([5]) is expressed as 



IV>o) = (cosf5 fe |0 fe ,0_ fc ) +isin0 fc |lfc, !_*)), 



(11) 



fe>0 



where the angles 9k satisfy the relation 2 ^ 
7 sinfc 



tan 26 



h + J cos fc 

7& 



h+ J(-l + fc 2 /2)' 



(12) 



FIG. 1: The phase diagram of transverse XY model. The 
vertical bold lines denote Ising transitions from the ferromag- 
netic phase to the paramagnetic phase , whereas the horizon- 
tal bold line stands for anisotropic phase transition between 
two ferromagnetic phases. Critical exponents associated with 
both the Ising and the anisotropic transitions are v = z = 1. 
The multicritical points (at J x — J y and h — ±1) are denoted 
by MCI and MC2, respectively; the exponents are v\ = 1/2, 
v<z = 1 and z mc — 2. 



Here the states \0k,0-k), and |lfc,l_fc) are respectively 
the empty and doubly occupied states of Jordan- Wigner 
fermions with momentum fc and — fc. 



A. Scaling analysis of the spectrum near the MCP 

As we discussed earlier multicritical points are char- 
acterized by two (or more) tuning parameters^. This 
affects the scaling form of various observablcs. The nat- 
ural tuning parameters in the present model arc Ai = 

. It is convenient to fix 
y = 1. Then dispersion 
can be rewritten in the 



Jx + Jy — h and A2 = J x ~ J\ 
the overall energy scale J = 
of the low energy excitations 
scaling form 



A fe a fc«-/(L|Air,L|A a |^) ; 



(13) 



where z mc = 2 is the dynamical exponent at the multi- 
critical point, vi — 1/2 is the exponent assosciated with 
the divergence of the correlation length when one ap- 
proaches the MCP and is the correlation length ex- 
ponent governing the cross-over from linear to quadratic 
dispersion at Ai = (see also Ref. [32|]). Physically the 
exponent v 2 governs the scaling of the healing length de- 
termining the decay of the longitudinal correlations^. 
The scaling function / is explicitly written as follows: 



f(x 1 ,x 2 ) = 




(14) 



Because the MCP in our case is the intersection of two 
critical lines one of which terminates at the MCP A2 = 
(the horizontal line on Fig. [TJ), the scaling function / 
should satisfy some constraints. In particular, at x 2 2> 



4 



Xx,l, the spectrum of excitation should reduce to that 
for the Ising critical point characterized by the exponent 
zn = 1: Afe ~ k Zit . This implies that 

f(x 1 ,x 2 )~x? {z -- z « ) (15) 

This is clearly the case since for v% = 1, z mc = 2, and 
2it = 1 we have f(xi,x 2 ) ~ a; 2 at large £ 2 . Likewise in 
the opposite limit of large x 2 — > the spectrum should 
be characteristic of the anisotropic transition occurring 
at finite value of x\ ~ 1 with a finite gap at X\ = 0. This 
implies that for x 2 *C x\ <C 1 



f(x\ 1 x 2 ) ~ const, 



(16) 



which is also the case. One can also infer the infor- 
mation about the scaling of the minimum gap near the 
anisotropic transition which occurs at x\ =i*«l, Con- 
sidering the limit x\ = x\ and x 2 <C 1 and noting the gap 
in the spectrum scales as X^^ from Eq. (fl~3]) we find 
that near the MCP the gap scales as 



A min (Ai,A 2 )~ lArl 1 



(17) 



Using the explicit exponents for our model we find that 
A m i n (Ai,A 2 ) ~ ^/|A 1 ||A 2 |. For the generic direction 
where Ai w A 2 = A this equation reduces to A m ; n ~ 
|A| 3 / 2 . Similar analysis can reveal relations between cor- 
relation length exponents and the asymptotics of the scal- 
ing function. 



B. Scaling of the geometric tensor 

Next we will analyze the scaling of the different compo- 
nents of the geometric tensor. We will first perform the 
analysis specific for the transverse field XY model and 
then discuss generalizations to generic MCPs. Follow- 
ing eq. (4) and the definitions of Ai, A 2 , we can rewrite 
Eq. (TEH) 



tan(20 fe ) 



A 2 fc 



Ai - fc 2 /2' 



(18) 



It is easy to see that the components of the geometric 
tensor in our model take the form 

Xap = i£(lfc|0;JO*><OA|0A,|l*>, (19) 

fc>0 

where | OA;) is the ground state wavefunction in the k state 
and \lk) is the excited state of the pair of fermions with 
momenta k, —k: 



|lfe) = sin0 fe |O fc ,O_ fc ) - icos0 fc |l fc , l- fc ). 



(20) 



Using explicit properties of the wavefunction, this geo- 
metric tensor can be rewritten as 

k>0 1 k>0 1 



<9tan20 fc <9tan20i. 



fe>0 



dX r 



(21) 



With help of Eq. (jT8|) we can compute all three indepen- 
dent components of this tensor and find their scaling: 



Xl2 



-y 

■XT, ^ 



Ajfc 2 , 
4L^ At ' 

fc>0 h 



X22 



-y 

4L ^ 



fc>0 



(Ax - fc 2 /2) 2 fc 2 



4L ^ 



A 2 fc 2 (Ax -fc 2 /2) 



k>0 



(22) 



Now, let us study each term of the tensor separately. 
First we focus on the MCP where Ai = A 2 = 0. In this 
case clearly the only nonzero component of the geometric 
tensor is x 22 ; 



X22 



k>0 



/V s 



L 
24' 



(23) 



Note that this scaling is exactly the same as across the 
Ising quantum critical points (vertical lines in Fig. [TJ . In 
this respect the geometric tensor evaluated right at the 
MCP docs not contain any additional information about 
the additional singularity in the system compared to the 
standard critical point. 

The situation changes dramatically if we consider the 
scaling of Xa,/3 close to the MCP. In particular, if we 
approach the MCP from a generic direction Ai rj A 2 with 
Ai > (to be specific we assume Ai = A 2 = A) we see that 
the other diagonal component of the geometric tensor xu 
becomes singular at the quasi-critical points occurring at 
2A = k* 2 = {2-Km) 2 L~ 2 , with m = 1,2,.., where only 
one momentum mode k* = V2X dominates the sum in 
Eq. j22]). There we find 



1 

k>0 



AJP 



1 (k*) 



1 



L {k*) 12 L{k*f 



(24) 



The lowest possible value of the dominant momentum for 
the quasi-critical point is k* = 2tt/L (corresponding to 
A = 2ir 2 /L 2 ). At this coupling we find the maximum 
value for the 1 1 component of the geometric tensor 



max 
All 



27T 



(25) 



The next local maximum for xu wm occur at the second 
quasi-critical point corresponding to k* = At:/L and it is 
easy to see that it will be down by a factor of 2 6 = 64 
compared to the first maximum. Once we go sufficiently 
far from the MCP all momenta will start contributing to 
the value of %n m Eq. ([22]) so that 



Xn(A, A) 



dk 



X 2 k 2 



2tt [(A-fc 2 /2) 2 + A 2 /c 2 ] 2 



7T 

4A 2 ' 



(26) 



It is interesting that if we naively connect the scalings 
([25]) with ([26]) we will get a mismatch. Indeed one can 
expect that when A becomes comparable to 1/L 2 the sys- 
tem size dependence in \n should disappear and the geo- 
metric tensor should become a function of A only with no 
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L dependence. This is indeed the case for usual critical 
points. However, here the situation is very different. At 
this value of A (corresponding to the shorter correlation 
length £i ~ 1/vA becoming comparable to the system 
size) the 11 component of the geometric tensor starts to 
oscillate as a function of A and these oscillations per- 
sist until the moment when the contribution of the sin- 
gle dominant momentum fc* corresponding to the quasi- 
critical point becomes comparable with that of all other 
momentum modes. Comparing Eqs. (|2"5|) with (f2l>| we see 
that this happens when the second, longer, correlation 
length £2 ~ 1/A reaches the system size. So we conclude 
that highly oscillatory behavior of xn occurs in the in- 
termediate values of the coupling l/L 1 ^ 1 < A < l/L 1 ^ 2 . 

It is worthwhile to mention that in case of the Hamilto- 
nian (8) with even boundary condition, the quantization 
of k changes to k = (2m + 1)tt/L in the P = sector—. 
This shifts the position of the quasi-critical points by a 
factor of 1/L 2 , but does not change the L 5 scaling of 
X™ ax . In the thermodynamic limit when the components 
of the geometric tensor arc determined by the contribu- 
tions from many momentum modes the boundary condi- 
tions become unimportant. 

Similar analysis can be performed on the other two 
components of the geometric tensor. In particular, one 
can check that the off-diagonal component also has a very 
sharp singularity near the quasi-critical point rapidly 
changing from large positive to large negative value both 
of which scale as L A . And finally the X22 component also 
has two (positive) maxima on the left and on the right of 
the quasi-critical point which scale as L 3 . Weaker singu- 
larities are observed in the vicinity of other quasi-critical 
points. As in the case of X n they persist until both cor- 
relation lengths become smaller than the system size. 

Next we will attempt to generalize our scaling results 
to the case of generic MCP which has two different cor- 
relation length exponents v\ and v 2 . We will always as- 
sume that v\ < i>2 (when v\ = v 2 the scaling results re- 
duce to those of a conventional QCP). For this purpose 
we will make the assumptions similar to those used in 
Refs. [9II33I]. Namely, we will still rely on the assumption 
that elementary excitations are created in pairs with op- 
posite momenta. In Refs. 2j]._27J it was shown that this 



assumption is not affecting the universal scaling results, 
which alternatively can be obtained through the analysis 
of the correlation functions, but significantly simplifies 
the derivation. Besides in gapped systems lowest energy 
excitations generically have quasi-particle nature^. Then 



X a/3 =$>|dxJfc)<fc|d A ,|0). 
k 



(27) 



Here denotes the state with two excited quasi- 
particlcs having momentum k and — k respectively. Note 
that 

1 (0|A„Er»a(r)|ft) 



(0\dx a \k) 



(28) 



Aa 2A fe 

For relevant or marginal operators the scaling dimension 
of the perturbation A Q J] r » a (r) is the same as the scal- 



ing dimension of energy (alternatively this perturbation 
is irrelevant). Therefore the scaling dimension of the ma- 
trix element (f2"5|) is the same as the scaling dimension of 
1/ A Q . This can be expressed introducing the dimension- 
less scaling function 



(0\dxjk) = ± /Q (i|A 1 r i ,L|A 2 r 2 ). 



(29) 



From this we can immediately infer the scaling of the 
geometric tensor in the thermodynamic limit when the 
summation over momenta can be substituted by integra- 
tion: 



W(A X ,A 2 ) = l-^—F a ^(\X 2 r-/\X 1 \^), 



where 



F a ,p{x)= I d qf a (q,q/x)fp(q,q/x). 



(30) 



(31) 



Note that we implicitly assumed that the integral over 
momenta converges in the infrared limit. Otherwise one 
can expect nonuniversal cutoff dependent corrections to 
the geometric tensor coming from regularization 2 ^. For 
our particular XY model the diagonal components of the 
scaling function F a ^ are given by the following integrals: 



F 22 (x) 



dq 



2 2 
q z x z 



2 71 



[q 2 x 2 + (l-q 2 /2) 2 } 
dq q 2 x 2 (l-q 2 /2) 2 
2^ [ g 2 x 2 + (1 _ g2/2 )2] 



(32) 
(33) 



If we would consider the curve lAil" 1 ~ IA2I" 2 then 
the argument of the scaling function F at p would be of 
the order of one and the scaling of the components of the 
geometric tensor will be given by the direct generalization 
of that for a QCP: 



(34) 



However, this setup does not describe a generic situa- 
tion. Because Ai and A2 are natural tuning parameters 
of the system it is more appropriate to analyze the scaling 
of the geometric tensor along the lines Ai = A cos 9 and 
A2 = A sin 9, where 9 is some fixed slope. In this case it is 
easy to see that the argument of the scaling function F a ^ 
goes to zero as A — > because v\ < v 2 . In turn this can 
affect the scaling of Xa,p- In our case this is indeed the 
case and at small x and cos 9 > the integral in Eq. ([3"T]) 
is determined by the interval q ~ v2 ± x, i.e. by the 
vicinity of the finite momentum critical points at A2 = 
and k = \/2Ai. The asymptotics of the functions F a ^ 
are then set by the properties of these finite momentum 
critical points for the anisotropic transition (see Fig. [1} . 
E.g. the asymptotics of F 22 is set from requiring that 
X22 ~ \X 2 \ du ^~ 2 . Combining this with Eq. ([30]) we find 
that at small x F 22 {x) ~ x, which is indeed the case. Sim- 
ilarly from the fact that the fidelity susceptibility along 
Ai direction diverges as 1 / 1 A2 1 at the anisotropic critical 
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point we find that Fu(x) ~ 1/x at small x. While these 
results are quite specific for the XY model the fact that 
they can be obtained from combining Eq. (|30| with the 
critical properties of the transition occurring at A2 — > 
and finite Ai is generic. 

For the path parametrized by Ai = Acos# and A2 = 
A sin 8 with cos 6 < there are no singularities occurring 
at A 2 — > because the critical line at A 2 = terminates 
at the MCP. In this case the scaling functions F a p(x) do 
not have any singularities at x — > and we recover the 
regular scaling (I3~H) . 

And finally at very small values of A such that the 
sum in Eq. (|27[) is dominated by the single momentum 
k* corresponding to the quasi-critical point we find that 
the scaling of the geometric tensor is determined by the 
scaling of the minimum gap at k*. Combining the scaling 
analysis with this observation we find that the maxima of 
Xa,p are determined by the scaling properties of the func- 
tions f a evaluated at k = k*. In general using Eq. (|29|) 
and for simplicity assuming d = 1 one can write the ex- 
pression for the geometric tensor in the form ([7]) with 



F a ,p = — x 

A a Afi 

2nn 2nn 
LlArh'LlAal- 



2im 



2nn 



L|Ai|^'L|A 2 |^ 



(35) 



For the system size smaller than the longer correlation 
length this sum can be dominated by a single most di- 
vergent term which determines the scaling of the geo- 
metric tensor. As we saw this is indeed the case for the 
transverse field XY model. 

Finally, we consider a path characterized by Ai ~ A 2 
for a general r > (sec appendix A). It can be shownil 
that quasi-critical points appear only for r < r c = 
2 where xn diverges as xn ~ L 4 / r+1 . In contrast, 
for r > 2, the scaling relation saturates to Xn ~ L 3 . 



C. Analysis of the geometric tensor through 
correlation functions. 



As we mentioned earlier the same scaling analysis can 
be repeated without making any assumptions about the 
nature of quasi-particles. Alternatively one can analyze 
the scaling of imaginary time correlation functions of the 
quench operator——: 



1 f°° 

X«,MAi, A 2 ) = j2 J rG Q/3 (Ai, A 2 ,r)dr, 



(36) 



where 

G a p{r) = (0\d Xa H(r)d Xfi H(0)\0) 

- (0\d Xa H(T)\0)(0\d Xl3 H(0)\0), (37) 

and 

d Xa H(r) = e HT d Xa He' Hr . (38) 



is the imaginary time Heisenbcrg representation of the 
operator d x H. Using our notations and translational in- 
variance of the Hamiltonian we find 



Xa,/3(Ai, A 2 ) 



(It 



\v a (r, T)vp(0,0)) c , 



(39) 



where the subscript "c" means that only the connected 
part of the correlation function should be taken into ac- 
count. One can check that all the scaling results of the 
previous section can be reproduced through the language 
of the correlation functions. For the XY model we give 
the details of the derivation of the correlation functions 
in the Appendix (B). In particular, we find that at the 
quasi-critical point 



Gn(r) ~cxp(L-V) 
which again yields 



^ max t b 
Xll ~ L 



(40) 



(41) 



D. Numerical analysis of the fidelity susceptibility 
along the fixed path. 

Having discussed the scaling of the geometric tensor, 
it is instructive to explicitly show the obtained depen- 
dencies. To be specific we will focus on the fidelity sus- 
ceptibility along a particular path given by h = 2J y 
fixed and J x being varied, which is equivalent to setting 
Ai = A 2 = A. The system crosses the MCP {J x = J y as 
shown in FigfT]). As we discussed earlier the fidelity sus- 
ceptibility is just the convolution of the geometric tensor 
with the derivatives d\ a /dX (see Eq. ((5J) so that in our 
case 



Xf(A) = xn(A, A) + x 22 (A, A) + 25R X i2(A, A) 



(42) 



Clearly the fidelity susceptibility is dominated by the 
most divergent component of the geometric tensor. 
Combining Eqs. (f22|) we find 



XF 



1 r-s (J y sin k + 2J y sin k) 2 



fe>0 



A 4 fc 



_L V — 



(43) 



fe>0 



Right at the MCP the scaling of xf ~ L coincides with 
that of X22 (see Eq. (f2"3"|) V As A increases the fidelity sus- 
ceptibility is becoming dominated by the xn component 
of the geometric tensor having very pronounced singular- 
ity Xf ~ L 5 at the first quasi-critical point A = 2ir 2 /L 2 . 
In Figs. [2] and [3] we show the numerically computed fi- 
delity susceptibility for the XY-model. In Fig. [5] we focus 
on a larger range of J x which contains both the Ising crit- 
ical point occurring at J x — — 3 (with the expected linear 
scaling of xf with the system size) and the MCP occur- 
ring at J x = 1. While Fig. [3] focuses on the vicinity of 
the MCP. We see that indeed xf has very pronounced 
non-monotonic behavior on the one side of the transition 
corresponding to Ai > 0. In Fig. Q] we show the scal- 
ing of the maxima with the system size and confirm L 5 
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asymptotics of the maximum of the fidelity susceptibil- 
ity together with linear in L scaling of the minima of 
Xf- As we explained earlier the oscillatory behavior is 
due to proximity to the critical line corresponding to the 
anisotropic transition (A2 = 0) where critical points oc- 
cur at finite momentum. Because v\ = 1/2 < V2 = 1 the 
generic direction Ai = A2 = A is effectively pushed to- 
wards the anisotropic transition, indeed | Ai | ^ ^> | A2 1 ^ 2 - 
In turn this results in emergence of the quasi-critical 
points characterized by anomalously small gap scaling as 
I A| 3 / 2 occurring when the allowed value of the momentum 
crosses the gapless point at the anisotropic transition line 
(Ai (2nn/L) 2 /2) (see Fig.(5)). 



le+06 
100000 




FIG. 2: The variation of \f with J x , as obtained numerically 
using eq. (21), for a system size of L — 100. We have fixed 
h = 2J y — 2. The peak near the Ising critical point at Jx — 



—3 scales as L, whereas, the maxima near the MCP ( J x 
shows a L° divergence. 



1) 




L 20° 



FIG. 4: The scaling behavior of Xf ax with the system size L at 
the first (xpi) an d second (XF2) maxima near the multicritical 
point. The dashed line corresponds to the fixed slope L . 
Inset: The first minima of xf, denoted by Xf'" i n Fig- El 
shows a linear dependence on the system size. 



0.16 




0.1 0.12 0.14 
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FIG. 5: Gap Ak with A for the momentum modes k = 0.1, 0.2 



le+06 




FIG. 3: Same as in Fig. [2] but near the multicritical point 
{Jx = 1)- The oscillatory behavior of the fidelity susceptibility 
is a signature of the presence of the quasi-critical points. Each 
of the maxima, denoted by XFj, j = 1, 2, .., scales as L 5 . 



J x J y • Th-6 minimum of energy gap occurs at different values 
of A for different k, thus giving the quasi-critical points. 



quench is closely related to the probability of exciting 
the system per unit volume, which in turn is related to 
fidclit y 21 i 22 . In the limit of large transverse field in the fi- 
nal Hamiltonian, when the ground state is all spins up (or 
down), the defect density becomes equivalent to the num- 
ber of wrongly oriented spins. In the same spirit, we next 
look at the defect density (n ex ) for a fast quench start- 
ing from the MCP. For a usual critical point the defect 
density following a fast quench of magnitude A starting 
exactly at the QCP is given b y 21 i 22 (see also Fig. [5]) 



|A| 



yd 



(44) 



This scaling can be obtained from e.g. the adiabatic per- 
turbation theory which states that (see Ref. p3j ) 



E. Defect density following the instantaneous 
quantum quench from the MCP. 

Recently it was understood that the density of de- 
fects (which is equivalent to the density of excited quasi- 
particles in the XY model) produced during the quantum 



E 



dX'{0\dy\k) 



< |A| / dA'xp(A'). (45) 



Then the scaling (|4"4"]) immediately follows from that of 
Xf- In our case the integral of the fidelity susceptibility 
diverges at small A since xf ~ Vl^l 2 therefore we need 
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FIG. 6: (color online) Plot of kink density n cx with the quench 
amplitude A for fast quench through critical (red line) and 
multicritical (black line) points, as obtained numerically, with 
a system size of L — 5000 spins, and h = 2J y = 2. Negative 
values of A correspond to quenching into the paramagnetic re- 
gion, while positive values indicate the ferromagnetic phase. 
Inset: Sudden jumps in kink density for the same system size 
is shown to occur very close to the multicritical point, when- 
ever A equals any quasi-critical value. No such jumps are 
observed when system is suddenly quenched into the param- 
agnetic phase. 
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FIG. 7: Variation of kink density n ex with system size L for 
fast quench through multicritical point, with A kept fixed at 
a very small value 0.0001, and h = 2J y — 2. 
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FIG. 8: Plot of excitation energy Q as a function of quench 
amplitude A for fast quench through the multicritical point 
into the ferromagnetic region, as obtained numerically, with 
a system size of L = 5000 spins, and h = 2J y — 2 fixed. 
Inset: Plot of excitation energy Q as a function of quench 
amplitude A for fast quench through the multicritical point 
into the paramagnetic region, 




FIG. 9: Plot of diagonal entropy density s as a function of 
quench amplitude A for fast quench through the multicritical 
point, as obtained numerically, with a system size of L = 
5000 spins, and h — 2J y = 2 fixed. Negative values of A 
correspond to quenching into the paramagnetic region, while 
positive values indicate the ferromagnetic phase. Inset: Plot 
of diagonal entropy density for small quench amplitude A. 
Oscillations are observed near the quasicritical points inisde 
the ferromagnetic region. 



to analyze the scaling of n ex separately. Let us observe 
that as in the case of xf the excitations arc dominated 
by the matrix element (0|9.\i|fc) which scales as l/|Ai|. 
The second observation we make is that the sum over mo- 
menta k is determined by the proximity to the anisotropic 
transition with shorter correlation length, i.e. k ~ vA. 
Combining these considerations we find that 

nex-IAI^ 1 . (46) 

The other simpler way to get this scaling is to note again 
that for the generic quench where A2 ~ Ai the spectrum is 
dominated by the direction with smaller value of v (since 
then lAil^ 1 | A2 1 ^ 2 - Thus the characteristic momentum 
above which transitions are suppressed is approximately 
equal to the inverse of the shorter correlation length 
k l/\\\ Ul . This characteristic momentum (in power 



of dimensionality) determines the scaling of the defect 
density. The situation qualitatively changes for a quench 
with Ai < 0. In this case the spectrum is dominated by 
the critical line at A 2 = and we expect n cx ~ |A| dl/2 . 
Numerical analysis of fast quench into the paramagnetic 
region (Ai < 0) shows n ex ~ |A| for large |A| (|A| ~ 0.1). 
(see Fig. 6). The anisotropy of the MCP reveals it- 
self, however, for the quenches with very small ampli- 
tude A ~ 1/i 2 , where only for the positive sign of Ai we 
observe step-like structure. These steps appear once a 
new quasi-critical point emerges when A ~ (27r?n) 2 /2L 2 , 
where m = 1,2,3....!/ — 1. As a result, the momentum 
mode k m = 2irm/L gets excited with high probability 
close to one leading to a sharp rise of the order 1/L in 
the defect density. In contrast, between any two consecu- 
tive quasi-critical points there is no significant change in 
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the total number of defects (see the inset in Fig. [6]) and 
the density does not appreciably change with A. Same 
oscillations can be revealed by fixing the quench ampli- 
tude A at a very small value and varying the system size 
(see Fig. [7]). There the non-monotonicity is related to the 
fact that the defect density decreases as 1/L with the sys- 
tem size between the jumps where the total number of 
excitations is approximately fixed. 

We can extend the above analysis to the case of 
heat densit y 21 ' 22 Q (or the excess energy above the new 
ground state) which, unlike be defined even for 

non-integrable systems not describable in terms of inde- 
pendent quasi-particles. For a fast quench through MCP 
into the ferromagnetic region, Q is expected to vary as 



\vi(d-\ 



A 



1.5 



(47) 



as is verified numerically (see Fig. [8]). Further, in con- 
trast to n exi Q does not show abrupt rise at the quasi- 
critical points very close to the MCP, even though its 
slope with respect to A increases slightly at these points. 
This is a consequence of additional suppression of the en- 
ergy by a factor of A(fc*) ~ (A;*) 3 . A similar quench into 
the paramagnetic region gives 



A 



(48) 



1 



which is expected if we use the exponents v — d = 
relevant in this case. 

We can also analyze the scaling of the diagonal entropy 
density 2 ^, 



f 

L 



^Pk log(p fc ) + (1 -p fc )log(l -p k ) 



Lk>0 



(49) 



where pk is the probability to be in excited state of mo- 
mentum ft. Unlike n cx , the diagonal entropy density is 
defined even if we break integrability: 



• s = -j^^Pn\og{p n ) 



(50) 



where p n is the probability to occupy n-th eigenstate. 
The diagonal entropy is equivalent to von Neumann en- 
tropy of the time-averaged density matrix. For any k 
mode, s k = - [pk log(pfc) + (1 - Pk) log(l - Pk)] is close 
to zero for pk ~ (i.e., for A < k 2 /2) or pk ~ 1 (i.e., 
for A > fc 2 /2), and assumes a finite value of order unity 
for pk ~ 1/2. This results in oscillatory behavior of en- 
tropy density inside the ferromagnetic near the quasi- 
critical points (see inset of Fig. |9|). On the other hand, 
coarse grained diagonal entropy density inside the ferro- 
magnetic region as well as the diagonal entropy density 
inside the paramagnetic region arc dominated by v = 1, 
which gives 



|A| 



ud 



|A|. 



(51) 



In our present model, numerical results suggest s ~ ,\ - 93 
for quench into the ferromagnetic phase, whereas s ~ 
A 1 ' 25 for a fast quench into the paramagnetic region. (See 
Fig. !|. 



On a related note, a recent studyii concerning slow 
non-linear quench across the MCP following a path 
parametrized by Ai = — | A2 1 r for a general r > has 
shown defect density to continuously vary with r as 
n ex ~ r~ r / 6 (1 < r < r c = 2) in the transverse field 
XY model, where A2 varies with time t as A2 = t/r, 
—00 < t < 00. A similar behavior of fidelity susceptibil- 
ity is also observed when we consider paths parametrized 
by Ai w A2, where Xii ax diverges as Xii ax ~ L 4 / r+1 for 
r < r c — 2 (sec appendix A). 



III. CONCLUSIONS. 

We have studied the scaling of the geometric tensor 
and the fidelity susceptibility near a quantum MCP us- 
ing the example of a spin-i XY chain in a transverse 
field. Our analysis shows that \f oscillates near the 
MCP, with the peaks, given by x™ ax ~ L 5 , occurring 
at the quasi-critical points. In turn these points are as- 
sociated with finite momentum anisotropic transition ter- 
minating at the MCP. The oscillations occur in the inter- 
mediate region of system sizes bounded between shorter 
and longer correlation lengths characterizing the vicin- 
ity of the MCP. Associated with these jumps we found a 
step-like behavior of generated defect density for sudden 
quenches of small amplitude starting from the MCP into 
the ferromagnetic region with a new jump occurring once 
a new quasi-critical point is crossed. In the thermody- 
namic limit these oscillations disappear and the defect 
density scales as |A| d,yi , where v\ is a smaller correla- 
tion length exponent characterizing the MCP. A simi- 
lar analysis of heat density gives a scaling of \ v ^ d+z ^) ; 
whereas diagonal entropy density shows oscillations near 
the quasi-critical points, finally scaling as s ~ \ V2d in the 
limit of large A. In contrast, behavior inside the param- 
agnetic region is always dominated by the Ising critical 
line, with v = z = 1. 

Using the adiabatic perturbation theory we generalized 
some of our findings to the situations of generic MCPs 
and showed how various scaling relations follow from the 
analysis of the scaling functions describing the transition 
matrix elements or equivalently the non-equal time corre- 
lation functions of the perturbation along the critical line 
characterized by a smaller correlation length exponent. 
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Appendix A: Appendix: scaling of the geometric 
tensor for a non-generic path 

In the main text we focused on the scaling of the geo- 
metric tensor along the generic direction Ai ~ A2. While 
this setup is indeed most natural, it is instructive to 
analyze the scaling along other directions, in particu- 
lar A 1 = cAj. For simplicity we choose c = 1. This 
scheme lets us approach the multicritical point following 
a generalized pathii, parametrized by the exponent r. In 
Rcfs. [35.36] it was shown that for a nonlinear quench in 
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time t near a usual critical point where A ~ t r the expo- 
nent vi associated with the amplitude Ai gets replaced 
by the product rv\ . A similar scenario occurs in our case. 
Indeed the scaling of gap can be written as 



A, 



k z 
k z 



|Ai| |A 2 | 



fcl/^l ' 

|A| |A| 



(Al) 



where in the last step we have taken A2 = A}^' = A. 
Similar transformations apply to the scaling of matrix 
elements. 

For our setup of the XY model in one dimension the 
spectrum gets modified to^£ 



(A2) 



which give*^, for r < r c = 2/2/^1 = 2, Afc ~ k 2 / r+1 
owing to the presence of quasi-critical points at 2A = 
(2-Km/L) 2 ^ r , m = 1,2,.... Whereas, for r > r c = 2, 
quasi-critical points are absent, and the scaling saturates 
to Afc ~ k 2 . Repeating the same analysis as in the main 
text we find that the maxima of say 11 component of the 
geometric tensor near the quasi-critical point now scale 
as 



- r *+i 



X11 = 



(A3) 



for r < 2. In the case of r = 1 we get back our previous 
result X11 ~ L 5 . However, for r > 2 the exponent in the 
L-dependence of Xa,P saturates, oscillation disappear and 
we get 



(A4) 



which describes the scaling of the maximum value of the 
fidelity susceptibility for the Ising transition: 



1 



A 2 



Xf 



L ^ k 2 (k 2 + A 2 ) 2 



(A5) 



occurring at A ~ k « 2n/L. This critical power of r c = 
v-2 1 V\ separates the regimes of dominance of anisotropic 
and Ising transitions on the scaling of the geometric ten- 
sor. 



Appendix B: Evaluation of correlation functions 

The Hamiltonian H at constant J\ + J2 = 1 can be 
written as 



- yEfe-W-^-^E^ ( B1 ) 



where Ai = 1 — h and A2 = 7. Therefore, 
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(B2) 



Using translation invariance of the system we have^i 
X«^(Ai,A 2 )= I drJ2r(V a (T,j)V (O,O)) c , (B3) 

where 



<y a (T,j)V(,(p,0)) e = (V a (T,j)Vp(0,0)) 

- (V a (T,j)){Vf,(0,0))> 



(B4) 



In the above equation we have defined A(t) = e Hr Ae~ Hr 
for any operator A. As we have seen above, X11 fixes the 
scaling of fidelity susceptibility at a quasi-critical point . 
Hence let us concentrate on evaluating (Vi (r, j)V\ (0, 0)) c . 
We have^i, 

(F 1 (i,r)y 1 (0,0)) c = K(r)^(0)) - (cr](r))(a z o (0)) 

1 



K(r)aS(0))-(m,) 



1 

4tt 
1 

4tt 



e ijk e -A k r 



Ajk „ — Al.t 

e e 



4tt 

cos2<f>(/c)dfc 
sin2$(fc)dfc 



jjk e - A k r dk 



— 7T 
2 



(B5) 



where $(fc) = l/2tan _1 [7sinfe/(/i + cosfc)], and (m 2 ) 
is the magnetization along the z direction. Taking into 
account that Afc and cos2$(fc) are even functions of k, 
and sin2<f>(A:) is an odd function of k, we can rewrite eq. 
(B5) as 



W(j,^i(0,0)) c 



1 

47T 



e AkT cos jkdk 



1 

47T 
1 

4tt 



n 2 



3 -A fc T 



cos jfc cos 2$(k)dk 



sin j k sm2$(k)dk 



(B6) 



Now, cosjfc ~ sinjfc ~ 1. Let us consider the first 
term 



1 

47T 



cos jkdk 



(B7) 



At a quasi-critical point, defined by ft, + cos k* = and 
7 = 1 — /i, the minimum value of Afc is Afc* ~ k* , 
whereas, Afc ~ k 2 for all the other momentum modes. 
Therefore 



where we have taken dk ~ 1/L and cos jfc ~ 1. 



(B8) 
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Let us now concentrate on the second term 

1 2 

„-A fc r. 



1 

4tt 



cos jk cos 2<£>(fc)<ifc 



(B9) 



At a quasi-critical point, 7 sin A;* /(/i + cosfc*) — ► 
—00, which makes cos [tan -1 sin k*/(h + cosfc*)]) = 
cos(— 7r/2) = 0. However, for the remaining k modes, 
cos (tan -1 [jsink/(h + cosfc)]) « cos(— fc) ss — 1 + fc 2 /2, 
and we have 



-E 

k^k' 



-k<T 



cos kj 



(BIO) 



since fc 2 » fc 3 . 

Now, let us consider the third term 



h = 



1 

4-7T 



T sin j'fc sin 2$(fc)dfc 



(Bll) 



Again, at a quasi-critical point, sin2$(fc*) ~ 1, whereas 
sin2$(fc) ~ fc for all k ^ fc*. Therefore we can write J3 
as 



(B12) 



Therefore combining the above terms we get, 
(F 1 (i,r)F 1 (0,0)) c ~-le- fc * 3T , 



and 



max 
All 



(B13) 



(B14) 



